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Abstract 



We present results from a quantum and semiclassical theoretical study of 
the p X y and p xx resistivities of a high mobility 2-D electron gas in the presence 
of a dilute random distribution of tubes with magnetic flux $ and radius R, 
for arbitrary values of kfR and F = e&/h. We report on novel Aharonov- 
Bohm type oscillations in p xy and p xx , related to degenerate quantum flux 
tube resonances, that satisfy the selection rule (kfR) 2 = AF(n + 1), with n an 
integer. We discuss possible experimental conditions where these oscillations 
may be observed. 

PACS numbers: 72.10Fk, 73.40.Qv, 73.40.-c 
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Transport in a two-dimensional electron gas (2DEG) in the presence of weak inhomoge- 
nous magnetic fields has recently been the subject of considerable interest, both experimen- 
tal and theoretical Qj. This situation has been achieved experimentally by gating the 
2DEG system with a type-II superconducting layer. Abrikosov vortices are then produced 
by applying an external magnetic field perpendicular to the plane of the layers. In the 
ballistic transport regime and for low fields, when the density of vortices is small, clear 
modifications to the Hall resistance in the quantum regime have been measured |§. Pre- 
vious theoretical studies of this problem were restricted to the asymptotic quantum regime 
kfR > 1 |] and classical limit kfR < 1 |, However, the experiments have covered the 
interesting intermediate kfR crossover regime, with F = 1/2. Here, kf is the Fermi wave 
vector and F = $/$ with, <3> — h/2 the flux quantum. 

In this paper we present a full solution to this problem for arbitrary kfR and F. Our 
results identify a series of novel quantum oscillations in the galvanomagnetic properties 
of the 2DEG, that appear to be within the reach of an experimental confirmation. These 
oscillations can be seen at intermediate ranges of kfR and F(> 1/2) values and are related to 
the Aharonov-Bohm (AB) effect. The intermediate ranges of kfR are experimentally already 
achievable (e.g. [T|,0]). At the end of the paper we discuss two experimental set-ups that have 
been suggested to produce larger values of F. Here we are interested in the experimental 
situation considered in || where the 2-D electrons move ballistically between the flux tubes 
and the dominant transport mechanism can be assumed to come from electrons scattering 
off individual flux tubes. Under these conditions, as a first approximation, we can apply the 
results of linear response theory in the Born approximation 0. These results are formally 
the same as those obtained with the Boltzmann equation ||. The weak nonlocal localization 
limit has already been considered experimentally and theoretically ||. 

There are three important physical contributions to the electronic transport properties of 
this system: (i) For finite R the Lorentz force that leads to an asymmetry in the scattering 
process, (ii) A diffractive force, relevant in the < kfR < 1 regime and first considered 
by Iordanskii 0, that also yields a transversal contribution to the transport, and (iii) the 
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standard AB contribution ||. The Iordanskii term in p xy , which is not taken into account in 
the differential cross section, is due to the scattering of electrons by finite radius flux tubes 
and has essentially the same origin as the AB effect U, for both are topological in nature 
and due to the long range properties of the vector potential. This means, as we see below, 
that the contribution from (ii) to the Hall resistance only depends on the value of F and 
not on the specific magnetic flux profile chosen in the analysis. 

The modification to the Hall resistivity, p xy , due to the inhomogenous field can be rep- 
resented by a Hall coefficient, a, which is defined by the expression || 

p xy = a(k f R,F) , (1) 

n e e 

where n e is the electron density and B the magnetic field. In the Born approximation of 
the Kubo formula the transport coefficients are expressed in terms of the scattering cross 
section /(</>), with the electronic scattering angle |6|,|]]. Explicit limiting values of a have 
been calculated in the extreme quantum a(kfR C 1) = s * n and semi-classical 

limits a(kfR 3> 1) = 1 [f^Q. We use these two results as constraints to be satisfied in our 
calculations. Previous studies were restricted to these two limits because of mathematical 
difficulties in the evaluation of the scattering amplitude f(<p) in the whole kfR and F 
ranges. These difficulties were identified by Khaetskii ||] and are essentially related to the 
singularity in the AB and Iordanskii scattering in the forward direction. In the AB case 
/ab(</> ~ 0) ~ l/sin(0/2), which would lead to an infinite a |5|,0]. 

Below we present the results of an explicit evaluation of /(</>) in the whole range of kfR 
and F values. More importantly, we use these results to calculate a and the magnetoresis- 
tance p xx in the extended parameter range. Since the calculational problems arise in the 
forward scattering region we consider the regularized scattering amplitude 

— J7T/4 OO 

/ e (0) = Z_ £ e^e-H e [e 2 ^-l], (2) 

m=—oo 



with e the regularization parameter, which is taken to zero at the end of the calculations. 
Here 5 m is the phase shift associated with the m partial wave and can be written as 



5 m = 6 m B — 8m- The 5^ B = |(|m| — \m + F\) accounts for the AB phase shift, and 
5 m = tan~ 1 (b m /a m ) for the remaining contribution to the scattering. The coefficients a m 
and b m are obtained from the asymptotic wave function solution to the Schrodinger equation 
ip(r — > 00,6*) ~ J2m=-ool a m Ju(kfr) + b m N u (kfr)} e md , which has the required form for 
incoming plane waves and outgoing circular waves. Here J v (kfr) and N u (kfr) are the Bessel 
and Neumann functions of order v — \m + F\ (see Ref. pi for more details). 



Including the contributions (i-iii) we can then write the Hall coefficient as, 

kr f 27T 1 

a = -f=]im sin0|/ e (0)| 2 C /0+— sin(FTr). (3) 

Lilt £->OJo nr 

The first term is the regularized Boltzmann contribution while we wrote the second topo- 
logical term following Iordanskii 0. An important property of this expression is that 
it fully reduces to the extreme quantum and semiclassical results mentioned above. In 
this expression, as long as e is finite, there is no singularity in f e {4>) and we can perform 
the integral. After evaluating this integral using Eq.(2), and from the general fact that 
[S m+ i — 5 m ) — > in the limit \m\ — ► oo, we get the finite result 



a 



sin[2(5 m+ i - 6 m )]. (4) 



This equation is one of the main results of this paper, for it provides an algorithm to calculate 
a for arbitrary values for the parameters kfR and F. The number of terms needed in the 
sum will depend on the parameter range considered. 

The corresponding linear response theory result for the magnetoresistance Ap xx is 



— - 1 = N F ii lim / "(1 - cos0) \hm 2 d<P (5) 



2 +oc 

= r N F £i sm 2 (6 m+1 -6 m ) (6) 

where r is the time between electronic flux tube collisions, within the Kubo-Born approxi- 
mation, and Np the concentration of magnetic flux tubes. When the magnetic field is zero 
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the transport scattering time is Tj = ii/vf, with Vf the Fermi velocity and £j the impurity 
scattering elastic mean free path, which is assumed to be much larger than 2R. 

We now proceed to present our results from the direct evaluation of Eqs.(4) and (6). 
Later in this paper we present a physical explanation of these results using a semi-classical 
analysis. In plotting these results we have used the typical experimental parameter values 
given in the captions. In Fig. 1(a) we show a as a function of kfR for different values of F. 
The curve F = 1/2 corresponds to an extended range of Fig. 3 in Ref. |§. We note that for 
values of F < 1, a is a monotonic function of kfR. Notice, however, that for F = 3/4, a 
can become negative for small values of kfR, which comes from our careful treatment of the 



extreme quantum region. For F > 2 we see clear oscillations in the a vs kfR curves [pL4 
For F = 10, for example, we can clearly identify sharp oscillations of a vs kfR, although 
their absolute value is smaller. The number of oscillations as a function of kfR is equal to 
the integer part of F, [F]. For F = 10, there are ten oscillations (five of them not shown 
occur for kfR > 15). Moreover, for small kfR there are narrow oscillations superimposed 
on the first few oscillations. In Fig. 1(b) we show the Hall resistivity as a function of F and 
kfR = 10. For small values of F we see the classical linear behavior of p xy up to a maximum 
value, after which it decreases as F increases. We note that the quantum curve, obtained 
using Eq.(4), decreases non monotonically as F increases and even becomes negative for 
values of F ~ 50. Finally, in Figs. 2 (a) and 2(b) we show the corresponding results for 
Ap xx (F)/ p xx (0) as a function of both kfR and F. In Fig. 2(a) for F < 1 we see that 
Ap xx (F)/p xx (0) is a monotonic decreasing function of kfR, while for larger values of F the 
magnetoresistance becomes an oscillatory function of kfR. In Fig. 2(b) we note the sharp 
resonances that occur exactly at the same values of the minima in p xy . 

We now provide a physical interpretation of these results in the semi-classical limit. To 



understand the semi-classical analysis, we start by discussing the classical problem [15 



We note that the energy, E = ^m*v'j, and the total (particle+field) angular momentum, 
J = m*Vfb — e$, are constants of the motion. Here b is the impact parameter and m* the 
electron's effective mass. The impact parameter is defined as b = y(t — > — oo), where y 



is along the direction of the current. The classical Hall coefficient is characterized by the 
important parameter (3 = uo c T = e§ /(2imi*VfR) = F/kfR. Different scattering events have 
different total angular momenta and different (3 parameters. When j3 < 1, the electron 
trajectories are only slightly affected by the magnetic flux tube. As (3 increases the Lorentz 
force becomes important until a critical /3 C , above which trapped orbits can exist. The 
particular quantitative value of (3 C depends upon the particular flux tube profile. For our 
flux-tube model (3 C = 1/2, and the trapped circular orbits have radius r Q = R/2/3. Both 
the quantum and classical scattering problems can be separated into angular and radial 
components. The radial component of the classical equation is, as usual, a one-dimensional 
problem with effective potential, V e ff(r,b) — [J + (3 J r r'B(r') dr'} 2 /2r 2 . Here the magnetic 
field of the flux tube is B(r) = B(r)z, and we have rescaled energies by m*vj, the angular 
momentum by 2nmVfR = hkfR, and distances by R. In these units J = b — (3 and the flux 
tube radius is equal to 1. In the inset of Fig. 1(b) we show curves for V e //(r, b) for different 
values of J. We see that electrons with different J's, or impact parameters, experience 
different effective potentials. As the total angular momentum decreases, V e ff(r,b) develops 
a potential barrier with height [J + (3} 2 /2, which decreases rapidly. We can show that for 
P > P c — 1/2 there is a range of J's for which there can be trapped circular orbits inside 
the flux tube from Ji = decreasing to J 2 = 1 — ft. As J decreases from J\ to J 2 , the 
center of the electronic circular orbit shifts from r 1 = to r 2 = 1 — r = 1 — ^. This range 
of possible total angular momenta is such that the circular orbit stays completely within 
the flux tube. Classically, these circular orbits can not be reached by a scattering process. 
However, quantum mechanically the scattering electron can tunnel through this potential 
barrier and form a quasi-bound state inside the flux tube for a finite time, and then escape 
again. In the classical calculation of the Hall resistivity and magnetoresistivity shown by 
dashed lines in the figures, we computed the classical differential cross-section which was 
used in Eqs.(l), (3) and (6), in place of \f{4>)\ 2 . 

In the semi-classical analysis we associate a classical circular orbit to each quasi-bound 
state. Using the standard Bohr-Sommerfeld quantization condition it is easy to show that the 
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quasi-bound states are degenerate and occur at quantized values of the energy, E n /Tiuj c = 
{kfR) 2 /AF = n + 1/2, with n an integer. This result is the analog of the Landau level 
condition in a homogeneous magnetic field. The factor n + 1/2 gives the total number of 
flux quanta enclosed by the circular orbit. Since the quanta of flux in the tube is equal to 
F, the quantum number n ranges from to [F] — 1. Moreover, the quantized total angular 
momentum, J m = mh, leads to a degeneracy of the n levels. This degeneracy is equal to the 
total number of quantized circular orbits which we can put inside the flux tube. From the 
range of classically allowed circular orbits mentioned above, we deduce that the allowed m 
values start at m\ = [{kfR) 2 / AF] = n and decrease down to m 2 = [kfR — F] + 5, with 5=1 
if [kfR — F] > and 5 = if [kfR — F] < 0. Therefore, we conclude that the degeneracy is 
equal to m x — m 2 + 1 = n + l + [F — kfR] — 5. In the figures, the arrows indicate the position 
of the principal quantum number n calculated using the selection rule {kfR) 2 = 4F(n + l/2). 
We observe that they are remarkably well aligned with some maxima and minima of p xx and 
p xy . Furthermore, we have numerically determined that each resonance in p xx and p xy occurs 
at a preferential angular momentum J m = mh for some m. We arrived at this conclusion by 
evaluating the time delay t^(kfR,F) = 2h(d5 m /dE) = (2R/vf)[dS m /d(kfR)} and we found 
that as a function of m, t 1 ^ becomes sharply peaked for one particular value of m each time 
the pair {kfR, F) corresponds to a resonance in the transport coefficients. In the figures we 
have indicated the values of m for the resolved resonances. The minima and maxima in p xy 
correspond to quasi-bound states due to the tunneling of the electron into the flux tube. 
Semi-classically, as m decreases the resonances correspond to rotationally asymmetric orbits 
leading to larger p xy , as observed in Fig. 1(b). Note that the number of resonances observed 
for a particular quasi-bound state level should be equal to the degeneracies of this level. 
However, near mi, the amplitude of the resonances is suppressed due to the exponentially 
small tunneling probability through the potential barrier (= [J + f3] 2 /2). For example note 
that p xy = at F = 50 in Fig. 1(b). On the other hand, for each quasi-bound state level, the 
observed resonances with the smallest m (m=3, 0, -7, -51 in Fig. 1(b)) correspond precisely 
to the value m 2 derived semi-classically. 
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We now consider the possible experimental conditions necessary to observe the galvano- 
magnetic oscillations described in this paper. The variation of kfR in the ranges of interest 
has already been achieved [1]J3|]. New techniques need to be developed to produce larger 
values of F inside the flux tubes. We discuss a couple of possibilities that have already been 
suggested to us. The general idea is to have the usual Hall bar shown schematically in the 
inset of Fig. 2 (a), with the inhomogenous magnetic field produced by a dilute distribution 
of magnetic flux tubes of strength F. One possible way to get larger values of F experi- 
mentally is by depositing randomly located submicron size superconducting dots or pillars 
on top of the 2DEG, in a manner similar to the way the dot and antidot systems have been 
fabricated Alternatively, one may drill randomly located sumbmicron holes in the 

superconducting layer by using electron beam lithography fi~3fl . In both cases, by following 



a magnetic field cooling technique the magnetic flux may be pinned inside the dots thus 
trapping a large bundle of flux quanta. As in the antidot systems we do not expect that the 
oscillations found here will be significantly affected by temperature or Coulomb effects, pro- 
vided the temperatures are sufficiently low and the charging energy effects are not significant 
for the superconducting pillars fabricated. 

In conclusion, we have presented a detailed analysis of the transport properties of a 2-D 
electron gas system in the presence of a dilute gas of randomly located magnetic flux tubes 
for arbitrary values of kfR and F. The main result from our analysis is the presence of novel 
AB-like oscillations in the galvanomagnetic properties of the system. These oscillations are 
explained in terms of the degenerate resonant levels, satisfying the selection rule (kfR) 2 = 
4F(n + |), due to the effective trapping potentials produced by the flux tubes. A more 



extensive presentation of the results described here will appear elsewhere ]T6 . 
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FIGURES 

FIG. 1. (a) Hall coefficient a as a function of kfR. Here R = 100 nm, the density of flux tubes 
is Np = 10 5 mm -2 and the electron concentration is n e = 3.98 x 10 10 cm" 2 , with <E> = irR? . Here 
n denotes the flux tube resonances and m the angular momentum degeneracies, (b) same as in (a) 
for the Hall resistivity p xy as a function of F. The dashed line corresponds to the classical results. 
The inset shows the effective potential V e ff f° r different values of the total angular momentum 
(from top to bottom J=1.25, 0.5, 0.2, 0, -0.25, -0.75) 

FIG. 2. Ap xx (F) I pxx{0) vs kfR (a) and vs F (b) for the same parameter values as in Fig.l, 
with mean free path £i = 2 pm. The inset shows the system considered in this paper formed by a 
Hall bar with a dilute random distribution of perpendicular magnetic flux tubes of strenght F&q. 
See text for further details. 
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